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Trapping and guiding surface plasmons in curved graphene landscapes 
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We demonstrate that graphene placed on top of structured substrates offers a novel approach for 
trapping and guiding surface plasmons. A monolayer graphene with a spatially varying curvature 
exhibits an effective trapping potential for graphene plasmons near curved areas such as bumps, 
humps and wells. We derive the governing equation for describing such localized channel plasmons 
guided by curved graphene and validate our theory by the first-principle numerical simulations. The 
proposed confinement mechanism enables plasmon guiding by the regions of maximal curvature, and 
it offers a versatile platform for manipulating light in planar landscapes. In addition, isolated defor¬ 
mations of graphene such as bumps are shown to support localized surface modes and resonances 
suggesting a new way to engineer plasmonic metasurfaces. 
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Introduction. Unprecedented optical properties make 
graphene a promising plasmonic material with great po¬ 
tential for practical applications ranging from optical 
routing [1-5] to nonlinear optics [6], optomechanics [7] 
and sensing [8] . In addition to relatively strong and tun¬ 
able electromagnetic response of this one-atom-thick ma¬ 
terial, it was shown that at infra-red frequencies graphene 
supports surface plasmon-polaritons (SPPs) with ex¬ 
treme confinement and propagation characteristics con¬ 
trollable by doping [9]. In particular, thanks to the strong 
dependence of its optical characteristics from electric bias 
graphene is on the path to be widely used for electro- 
optical modulation [1, 10]. Possibility of such electro¬ 
static doping enables control over the electrons density 
at the Fermi level which defines high frequency optical re¬ 
sponse of graphene [11-15]. With tremendous success in 
graphene electronics witnessed recently, these unique ca¬ 
pabilities enable integration of plasmonic and electronic 
devices on the same substrate [16-20] and even a design 
of tunable cloaking devices [21]. 

Presently, two of the most common approaches to uti¬ 
lize unique plasmonic response of graphene for guid¬ 
ing applications rely on (i) inhomogeneous gating of 
graphene [22, 23] affecting its local plasmonic response, 
and (ii) subwavelength-scale patterning of graphene layer 
giving rise to localized plasmonic resonances [5] . 

In this Letter, we propose an alternative approach to 
engineer plasmonic guiding properties of graphene by 
employing its property to confine and scatter light by 
curved regions such as extended one-dimensional humps 
and local two-dimensional deformations. Technologically, 
such curved graphene profiles, shown schematically in 
Figs. l(a,b), can be fabricated by either transferring 
graphene onto curved dielectric substrates or grown di¬ 
rectly by laser ablation on structured surface of SiC 
demonstrated recently [24]. In particular, we demon¬ 


strate plasmon guiding by a bump of monolayer graphene 
which effectively operates as a plasmonic channel waveg¬ 
uide [25-33]. In addition, we show that a local isolated 
deformation of graphene in the form of a bump supports 
local surface plasmonic resonances. 

It has been already established that a curved graphene 
surface can play a role of an effective potential for elec¬ 
trons and phonons, and it may support spatially local¬ 
ized modes [34]. Similarly, we may expect that graphene 
plasmons can be modified by effective geometric poten¬ 
tials, as this was shown earlier for curved metal-dielectric 
interfaces [35]. We notice that the study of the propa¬ 
gation of light in curved space has attracted some spe¬ 
cial attention due to novel opportunities to guide and 
focus light [36-38]. Here, we extend this knowledge to 



Figure 1. (Color online) Schematic view of convex curved 
substrates covered by graphene and associated global Carte¬ 
sian and local curvilinear coordinate systems, (a) and (b) 
correspond to the cases of one- and two-dimensional curved 
surfaces, respectively. 
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the case of plasmons in curved graphene and demon¬ 
strate novel approaches to engineer plasmonic response 
of curved graphene. 

Analytical approach. As the first step, we explore 
the capability of trapping light with the use of an an¬ 
alytical approach assuming graphene surface of small 
one-dimensional curvature. Specifically, we consider a 
graphene monolayer placed on a curved interface between 
two dielectric media with permittivities £_ and £+, as 
shown in Fig. 1(a). Such a graphene-coated surface forms 
a channel waveguide for surface plasmons, thus enabling 
a two-dimensional (2D) spatial localization of light, al¬ 
ternative to waveguides created by graphene conductiv¬ 
ity modulation [22, 23]. We base our analytical consid¬ 
eration on Maxwell’s equations solved in the adiabatic 
approximation within the small-wavelength limit. 

We assume that an inhomogeneity in the graphene pro¬ 
file is described by a smooth function x = /(z), and we 
look for modes propagating along the y axis. For con¬ 
venience, the curvilinear orthogonal coordinate system 
{£, 2 /, s} bound to the reference curve x = f(z) is used in 
place of the Cartesian coordinate system {x,y,z}. Ac¬ 
cepting harmonic exp(— iujt) time-dependence, Maxwell’s 
equations can be written as follows 

( V x E = iko H , 

< 47r (!) 

V x H = -ik 0 £(g)V + —6(£)(t(u>) E r , 

1 c 

where fco = cj/c is the wavenumber in free space, cf(uj) = 

+a( R \oj) is the linear frequency-dependent sur¬ 
face conductivity of graphene, the function e(g) describes 
the dielectric permittivity distribution, 


<o = 



£<o, 

£> 0 . 


( 2 ) 


The Dirac delta-function 6(g) in Eq. (1) indicates that 
the graphene layer is placed at g = 0, and the subscript 
r refers to the field component tangential to the surface. 
Equation (1) leads to 

Ajr 

-VxVxH+*ge(OH = *A;oVe(OxE T - <t(w)Vx< 5(£)E t , 

(3) 

where E r = -——— (V x H) is continuous at g = 0, and 
k 0 £(g) 

the vector operations can be calculated using the Lame 
coefficients given as = 1, h y = 1, h s = 1 — >c(s)g = 

1 — g/R, where x(s) ~ f"(z ) is a signed curvature of the 
reference cylinder line x = f(z), R = 1/ k(s) is the local 
radius of curvature. 

Following the approach employed earlier in Refs. [39- 
41] , we derive the equation for the slowly varying plasmon 
amplitude, employing the asymptotic description often 
used in the paraxial optics and physics of optical soli- 
tons [42] . To develop the consistent perturbation theory, 
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Figure 2. (Color online) Top panels in (a,b): Effective po¬ 
tential (black solid curves) and transverse profiles a(z) (a.u.) 
of the trapped modes (blue dotted curves) calculated for the 
Gaussian profiles f(z ) = ±/o exp (— z 2 /w 2 ), fo = 200 nm, 
w = 506 nm, and dielectric permittivities e+ = 1, e_ = 2 
at the frequency cJo/27r = 16 THz. The red dashed line lev¬ 
els the relative correction A/c sp //c sp (cl;o) to the wavenumber 
fep(^o) = 22.4/co corresponding to the fundamental channel 
mode. Bottom panels in (a,b): Normalized magnetic field dis¬ 
tribution | H | /1H | ma x of the modes shown in top pannels and 
obtained with the use of first-principle FEM solver. 


we introduce a small parameter, 


fx 2 = max 


1 

ko\R\ 


cr( J ) 


assuming the losses in graphene and spatial inhomogene¬ 
ity to be small. Accordingly, the solution of Eq. (3) is 
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sought in the form 


where the damping coefficient 7 is expressed as 


H* = 


Hy = 


A(/U 2 y,A<s)MO + /U 2 H (2) (A« 2 y,A<s,e + 


0 ik S py 


A(s,y)h(0 + H^(s,y,0 + - 


Jk sp y 


MH ( 1 ) (M 2 y>MS,C) + m 3 H ( 3 ) (m 2 Y>MS,C) + ••• 


0 ik sp y 


H^(s,y,0+H^(s,y^) + ... 


rjikspU 


( 4 ) 

where A(s,y) is slowly varying amplitude of the plas- 
monic mode propagating along the curved channel and 
k sp (oo) satisfies the dispersion relation 


£+ £- .47r . x 4 tt m , N 

-t-= —7 — cr(u;) = —, 

K+ K- 00 00 


1/2 


( 5 ) 


where /^+,- = (fcg p — &o£ +j _) 

Substituting the expression (4) into Eq. (3) in the zero 
order in /x returns the dispersion relation (5) and trans¬ 
verse profile for the plasmon on a locally flat graphene, 


HO = -ik 0 < 


^± e -*+€, 


- e 

K- 




e >0, 
£ < 0 • 


(6) 


The correction of the first order in /x is determined 
from V • H = 0 as 

In the order of /j 2 the correction H {2) at £ 7 ^ 0 satisfies 
the equation 


d 2 H ( 2 ) 

. V^vs, -sp 


+ (fc 0 2 e(C - fc s 2 p )il (2) = F 
dA d 2 A 


F — — ( 2 zfc sp —-K 


3s 2 


/ j.\ , x „ dh 

MO - x(s)^l— 


( 7 ) 


Matching solutions of Eq. (7) at the interface £ = 0 by 
using the field boundary conditions 


7 2) (£ = +o) = 7 2) (e = -o) = 7 2) , 

47T , 
c 

where 


H [ 7c = +0) -17 (2) (£ = -0) = -^ (a ( 7wM+^ (/ V)7 2) )» 


4 2 7 = ±0) = .4 = ±0) - , 






finally we find that the amplitude A of the plasmon sat¬ 
isfies equation 


2«fc S p ( -7^- + 7 A) + 


cPA 
ds 2 


-5(«M = 0 , 


(8) 


CJ 

Q T Q J rusp 


and the curvature-dependent function g(s), 


g(s) = k(s) 





is an effective geometric potential for plasmons [35] . 

Remarkably, the result of Ref. [35] for plasmons at 
a curved metal-dielectric interface in the absence of 
graphene can be retrieved from Eq. ( 8 ) by setting a(oo) = 
0 , and assuming that the substrate possesses a plasma- 
like frequency-dependent dielectric permittivity In 
this case £_//€_ = — £ + /k + , that leads to g(s) = 

-R- 1 = R- 1 ^- ( ---) ' . Besides, a 

K+ — K- ko \ S+ + S- J 

particular case characterized by a constant curvature re¬ 
covers the problem of TM plasmon propagation along 
cylindrical graphene-coated nanowires [43]. 

It is also important to note that, as follows from 
Eq. ( 8 ), the symmetric dielectric environment, £_ = 
£ + = e, suggests g(s) = 0 , and implies that no local¬ 
ization to the curved region is possible in the case of 
ID curvature, while the loss coefficient is found to be 



Figure 3. (Color online) Channel waveguide modes found 
for the value of graphene conductivity a — 0.00035x S using 
FEM method for = 1, &- — 2: shown is the surface charge 
density distribution p s on the curved graphene surface. The 
frequencies of the eigenmodes plotted in (a) and (b) are 21.1 
THz and 23.9 THz, respectively. The geometrical parameters 
of the curved region are: the height (depth) of the channel is 
300 nm and the width of the upper (lower) elliptical segment 
is 200 nm, and the radius of curvature of the two lower (upper) 
circular segments is 140 nm, shown y -length is 1 /xm. 
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Figure 4. (Color online) Geometry and localized modes of (a) 
circular channel (depth 100 nm, side curvature radius 33 nm) 
and (b) hemispherical bump (height 200 nm, side curvature 
radius 66 nm) found with the use of FEM method. Upper pan¬ 
els in (a) and (b) show profiles of the curved regions. £+ = 1, 
£- — 2. Four lower panels in (a) and (b) depict the surface 
charge density distributions corresponding to the modes with 
increasing number of nodes in angular and radial directions. 
The eigenfrequencies found for the value of graphene conduc¬ 
tivity a = 0.00035z S (from top to bottom) (a) 11.1 THz, 
14.4 THz, 15.2 THz, 21.8 THz; (b) 10.8 THz, 14.5 THz, 17.2 
THz, 22.8 THz. 


7 = 


k S pid 


t( r \u) — , which is in agreement with that 


calculated earlier in Refs. [40, 41], where the equation for 
the correction H^ is written with incorporated bound¬ 
ary conditions using the formalism of S functions. 

As an example, we consider the interface profile of the 
Gaussian form. This results in inhomogeneous curvature 


which at some parts becomes a potential well capable 
of confining a mode. In our calculations, graphene is 
described by the surface conductivity written in terms 
of the Drude model cr(u) = (ie 2 /7rh 2 )£F(uj + ^RntnJ -1 , 
where the Fermi energy Ep — 0.307 eV, Ti ntra = 0.1 ps is 
the relaxation time. This case of doped graphene satis¬ 
fying hu < S F , allows neglecting the effects of interband 
transitions and we neglect temperature effects [12]. 

The potential corresponding to the Gaussian chan¬ 
nel operates as a 2D plasmonic waveguide with trapped 
modes being eigenstates of the stationary Schrodinger- 
like Eq. (8), shown in Fig. 2. The first example consid¬ 
ered here corresponds to the substrate dielectric constant 
exceeding that of the superstrate. Only one bound state 
is found for the given set of parameters with its amplitude 
spatial dependence A(z, y ) = a(z ) exp (iAk sp y). The po¬ 
sition of eigenvalue A k sp with respect to the trapping 
potential is indicated in Fig. 2 (a) (top panel, red dashed 
line). The state appears to be less localized to the curved 
region with the fields primarily located near the side re¬ 
gions of maximal positive curvature. The second case 
studied has inverted profile and also supports one local¬ 
ized state illustrated in Fig. 2(b) top panel. One can 
see that for this case the field is mainly localized at the 
central dip (notch) having the maximal positive curva¬ 
ture and exhibits quite strong localization in the in-plane 
transverse z direction. 

Numerical approach. To validate the developed ana¬ 
lytical theory, first principle numerical simulations are 
performed with the use of commercial full-vector FEM 
solver COMSOL Multiphysics. The graphene is modeled 
by the surface current tangential to the curved interface 
between the substrate and the superstrate j = <jE t . The 
numerical results, shown in Fig. 2 (a,b) bottom panel are 
found to be in excellent agreement with the results of our 
analytical theory. 

While the perturbative nature of the analytical model 
limits its applicability to the case of small curvatures, the 
possibility to confine light with the curvature extends be¬ 
yond this case. Moreover, one can expect that the trans¬ 
verse in-plane confinement of the channel modes can be 
further improved by increasing curvature. The particu¬ 
lar cases of guiding by the strongly curved non-analytic 
graphene profiles are given in Figs. 3(a,b). Similar to 
the case of small curvatures, the mode corresponding to 
graphene curved towards the high index material shows 
more localized mode at the center of the curved region. 

The proposed concept of trapping light in curved 
graphene can be further extended to the case of 2D pro¬ 
files. Here we consider the cases of generalization of the 
ID profile closed to form circular channel and hemispher¬ 
ical bump in 2D geometry. The resultant structures pos¬ 
sess cylindrical symmetry and the corresponding trapped 
modes are expected to be quantized with respect to the 
azimuthal direction. Indeed, numerical results obtained 
in COMSOL by solving full 3D FEM problem for both 
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cases shown in Figs. 4(a,b) reveal a variety of modes with 
increasing number of nodes as the frequency of the mode 
gradually increases. As opposed to the channel modes, 
and similar to the case of localized surface-plasmons in 
metallic nanoparticles, these modes are leaky, and there¬ 
fore they have a finite radiative lifetime. Nonetheless, 
just as in the case of plasmonic oligomers and periodi¬ 
cally arranged metasurfaces, we expect that by building 
arrays of graphene bumps and dents one can engineer 
modes of high quality factors limited only by internal 
Ohmic loss in graphene. 

Conclusions. We have suggested a novel approach to 
confine light in curved graphene landscapes. We have 
shown that curved graphene allows guiding surface plas- 
mons trapped to the curved regions. With the technique 
to either directly grow or transfer graphene on top of 
structured substrates fully developed, the demonstrated 
ability to confine and manipulate light with the curva¬ 
ture envisions photonic networks of interconnected chan¬ 
nel waveguides, which can serve as a versatile platform 
for on chip graphene plasmonics integration. 
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